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ABSTRACT 

We present the results of Keck/NIRSPEC spectroscopic observations of three Lya emitting galaxies (LAEs) 
at z 2.3 discovered with the HETDEX pilot survey. We detect Ha, [O III], and H/3 emission from two galaxies 
at z = 2.29 and 2.49, designated HPS 194 and HPS256, respectively, representing the first detection of multiple 
rest-frame optical emission lines in galaxies at high-redshift selected on the basis of their Lya emission. We 
find that the redshifts of the Lya emission from these galaxies are offset redward of the systemic redshifts 
(derived from the Ha and [Olll] emission) by Av = 162 ± 37 (photometric) ± 42 (systematic) km s"' for 
HPS194, and Av = 36 ± 35 ± 18 km s"' for HPS256. An interpretation for HPS194 is that a large-scale 
outflow may be occurring in its interstellar medium. This outflow is likely powered by star-formation activity, 
as examining emission line ratios implies that neither LAE hosts an active galactic nucleus. Using the upper 
limits on the [N ll] emission we place meaningful constraints on the gas-phase metallicities in these two LAEs 
of Z < 0.17 and < 0.28 Zq (Icr). Measuring the stellar masses of these objects via spectral energy distribution 
fitting (~ lO'" and 6 x 10** M0, respectively), we study the nature of LAEs in a mass-metallicity plane. At 
least one of these two LAEs appears to be more metal poor than continuum-selected star-forming galaxies 
at the same redshift and stellar mass, implying that objects exhibiting Lya emission may be systematically 
less chemically enriched than the general galaxy population. We use the spectral energy distributions of these 
two galaxies to show that neglecting the contribution of the measured emission line fluxes when fitting stellar 
population models to the observed photometry can result in overestimates of the population age by orders of 
magnitude, and the stellar mass by a factor of ^ 2. This effect is particularly important at z > 7, where similarly 
strong emission lines may masquerade in the photometry as a 4000 A break. 
Subject headings: galaxies: evolution 



1. INTRODUCTION 

Star-forming galaxies at high-redshift are one of the most 
useful probes of the distant universe. By studying their physi- 
cal properties, we learn about the stellar mass, dust and chem- 
ical evolution during the first few Gyr after the Big Bang. 
The most efficient method of learning these properties is by 
comparing the observed colors of galaxies (their spectral en- 
ergy distributions; SEDs) to model stellar populations, finding 
which combination of metallicity, star formation history, ex- 
tinction, age and stellar mass best m atches the observed galax- 
ies (e .g., Papovich et al. 2001; Shap lev et al.l2001l:lGiavaliscol 
I2OO2I: iBruzual & Chariot 2003). This technique requires only 
broadband photometry, thus many galaxies can be studied 
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with a single dataset. 

SED fitting has been used to study thousands of galax- 
ies at high redshift, the majority of wh ich are either color 
select ed (Lyman break galaxies; LBGs; ISteidel & HamiltonI 
1 1 9931) . or selected on the b asis of a bright Lya e mission line 
(Lyg emitters; LAEs; e.g.. lCowie'& Hu 1998; R hoads etall 
120001) . Over 3 < z < 6, LBGs have been found to be 
highly star-forming galaxies (star-formation rate ^ lO's of 
Mq yr"'), with significant dust attenuation (Ay ^ 0.2 - 1.0 
mag), stellar population ages of lOO's of Myr and stellar 
masses of ^ 1 0^° - 10" Mp, (e.g ISawicki & Yed [19981: 
Papovich et all l200Tt [Shaplevetal.1 1200 1'. '2005^. 'Yan etall 
2005L 120061: lEvles et all I2005L l2007n 'Fontana et al. 200i 



lOverzier et al.1 120091; 



Reddv et al. '2006t IHuang etaP l20a 
Stark et al. 2009). 

Over the same redshift range, LAEs appear to be less 
evolved in the same characteristics. Their star-formation 
rates are typically more modest (< 10 Mq), they typi 
cally have ages < 100 Myr, and similarly masses of 
10'' Mq (e.g.,|Gawiseretal. 2006, 2007; Pirzkal et al. '20Q| 
Finkelstein et al. 2007, 2008, 2009d; Gronwall et al. 2007^ 
Lai et al..2007.,2008;,Pentericci et al.„2009;,Ono et al...2010bl 



< 



I Yuma et aTT 2010I) ." Recent results imply that at z > 7 
most galaxies have masses and ages similar to LAEs at 
z < 6, and more evolved, LBG-like galaxie s are rare ( e.g., 
Ouchi et al. 2009; Finkel stein et all I20l6bl: iGonzdlez etall 
i2010; Ono et al., 2010a) . While the physical properties of 
LBGs evolve from z = 6 - 3, the characteristics of LAEs 
appear unchanging, implying that they may represent newly 
forming galaxies at each redshift. 
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One of the more intriguing results to arise from high- 
redshift LAE stellar population analyses is t he result that 
many LAEs appear to have dust extinction (e.g.,'Pirzkal et alj 
2007; Finkelstein et al. 2008, 2009d; Pentericci et al. 2009; 
Lai et al .1120071, see also Blanc et al. 2010). This is surprising, 
as dust can efficiently attenuate Lya photons, as they reso- 
nantly scatter with neutral hydrogen, and thus can have long 
path lengths prior to escaping a galaxy. A number of sce- 
narios have been discussed in the literature explaining how 
Lya can escape through a dusty medium, including shift- 
ing out of resonance due to outflows (e.g., Ahn et al. 2003; 
Verhamme et all feOOSt iDiikstra & Loebi BOOl IZheng et alj 
2010) . as well as scatte ring off the surfaces of dusty H i 
clouds in a clumpy IS M (lNeufeldlfT99Tl; [Hansen & Ohll2006l; 
[Finkelstein et all 120081 12009dl) . These scenarios, or some- 
thing similar, are likely in play, as observations of Lya in local 
star-forming galaxies have shown that Lya can be observed in 
emission even though dust is present (e.g.. Atek et al. 2008.) . 
and both low- and high-redshift observat ions have found evi- 
dence that the Lya emission is extended jNilsson et alj|2009t 
lOstUn et al . 2009; Fi nkelstein et"ai]l2010cl) . 

Nonetheless, the level of uncertainty on the dust content 
in high-redshift LAEs is large, as the presence of dust has 
been deduced by fitting the SEDs of LAEs to stellar popu- 
lation synthesis models. In addition, these models place no 
constraints on the metallicities, as they suffer a number of de- 
generacies as a red spectral slope can be explained by dust, 
an older stellar population, or a higher metallicity. In order 
to learn about these crucial evolutionary properties in more 
detail, rest-frame optical spectroscopy is needed. However, 
at such high redshifts, these lines are shifted into the near- 
infrar ed, an d thus a re observationally challenging. 

Erb et al.l (120061) performed a large near-infrared spectro- 
scopic survey of star-forming galaxies at z 2.3, selected via 
the BX method to be similar intrinsically to LBGs at z ^ 3 
(ISteidel et al.ll200l . using the K eck telescope with th e near- 
infrared spectrograph NIRSPEC dMcLean et alj|1998l) . They 
primarily observed in the K-band, which allowed observations 
of Ha and [N ll] emission at this redshift (though they ob- 
tained H-band spectroscopy, and thus H/? and [O III] measure- 
ments for a few objects). From the ratio of these lines, they 
derived the gas-phase metallicity of these objects. Combining 
this information with stellar mass estimates, they examined 
the mass-metallicity relationship at z ^ 2.3, finding that star- 
forming galaxies at that redshift followed a similar sequence 
of higher meta llicity with higher mass as at very low red- 
shift (T remonti et al.ll2004 ). though the sequence was shifted 
down by ~ 0.3 dex in metallicity, showing that high-redshift 
star-forming galaxies are less chemically evolved than lo- 
cal star-forming galaxies. A number of other studies have 
taken advantage of magnification due to gravitational lensing 
to investig ate high-redshi ft LBG-like galax ies in detail (e.g., 
iTeplitz et al . 2000; Finke lstein et al.ll2009cl) . However, until 
now, the physical properties of galaxies selected on the basis 
of Lya emission have not been spectroscopically probed in 
detail. 

Observations of rest-frame optical emission lines can also 
allow a determination of the systemic redshift of a galaxy. 
Comparing this redshift to that of the Lya emission can al- 
low one to diagnose the kinematics of the inters tellar medium. 
Using a composite of ^ 800 LBG UV spectra, iShaplev et alj 
(|2003) found that Lya was redshifted by 360 km s"^ with re- 
spect to the systemic redshift, and the ISM absorption lines 



were blueshifted by 150 ± 30 km s"'. ISteidel et all ( 120101) 
studied the ISM kinematics of a sample of 89 galaxies at 
z ^ 2.3, finding that Lya and the ISM absorption lines were 
shifted by +445 and -164 km s~', respectively, when com- 
pared to the systemic redshift as measured by the Ha emis- 
sion. These velocity differences are thought to be due to the 
presence of global outflows in the ISMs of these galaxies. In 
such a scenario, the non-resonant nebular lines are observed 
at the systemic redshift, as they originate from the H II regions 
within the galaxy. The ISM absorption lines are blueshifted 
with respect to the systemic redshift, as the absorption fea- 
tures are created when the outflowing gas on the near side 
(with respect to the observer) absorbs stellar light. Lastly, Lya 
preferentially escapes after it has back-scattered off of the far 
side of the expanding ISM, shifting the Lya-line center out 
of resonance with H I, enabling it to traverse back across the 
galaxy, and through the near side of the ISM. 

Until recently, LAE samples were not available at 2 < z < 
2.7, which allows [Olll] and Ha observations with H- and 
K-band spectroscopy, respectively. LAEs are typically much 
fainter than LBGs, thus their emission line fluxes are fainter 
as well, resulting in more time-demanding observations. Re- 
cently, McLinden et al. (2010) published the first rest-frame 
optical emission line detections from LAEs, detecting [O III] 
emission from two LAEs at z 3. 1 . They were able to probe 
the kinematics of the ISM by comparing the redshifts of these 
lines to that of Lya, finding that both galaxies exhibited large 
scale outflows of 343 ± 18 and 126 ± 17 km s"'. Addition- 
ally, iHaves et"an (120101) probed the Lya escape fraction in 
high-redshift star-forming galaxies by measuring the Lya and 
Ha fluxes with optical and near-infrared narrowband filters, 
though only a few objects were simultaneously detected in 
both lines. However, probing further into the physical char- 
acteristics of such galaxies requires more than a single rest- 
frame optical line. We present the first detections of multi- 
ple rest-frame optical emission lines from galaxies at high- 
redshift selected on the basis of their Lya emission, which we 
use to study in detail their physical properties. In §2, we dis- 
cuss the selection of our LAE sample, and our near-infrared 
spectroscopic observations with Keck/NIRSPEC. In §3, we 
present our emission line measurements. In §4 we measure 
the systemic redshift and search for the existence of outflows 
in the ISM, in §5 we probe for AGN signatures, and in §6 we 
place constraints on the dust extinction. In §7, we constrain 
the gas-phase metallicities of our LAEs, including studying 
them on a mass-metallicity plane. Lastly, in §8 we examine 
how the measured emission line fluxes can affect the SED fit- 
ting results. We assume Hq = 70 km s~' Mpc~\ fini = 0.3 
and iljn = 0.7, and unless otherwise specified we use the AB 
magnitude system throughout (Oke & Gunn 1983) . 

2. DATA 

2.1. Sample Selection 

The Hobby Eberly Telescope Dark Energy Experiment 
(HETDEX) is a blind integral field spectrograph (IFS) search 
for LAEs at z = 1 .9 - 3.5, with a primary science goal of prob- 
ing dark energy via the power spectrum of 0.8 million LAEs. 
The survey will begin in 2012, using the Visible Integral-field 
Replicable Unit Spectrograph (VIR US) instrumen t, which is 
composed of 150 individual IFSs ( iHill et alj|2008b). Cur- 
rendy, a single prototype IFS (VIRUS-P; H iU et alj (l2008a)) 
is mounted on the 2.7m Harlan J. Smith Telescope at the Mc- 
Donald Observatory. Over four years, VIRUS-P has been 
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used for 111 nights as part of a HETDEX pilot survey, discov- 
ering 103 LAEs at z - 2-3 jAdams etal.ll2010t iBlanc et alj 
120101 hereafter Paper I and Paper II, respectively). The pilot 
survey probes a large volume relative to narrowband imaging 
searches, thus the discovered LAEs are bright, and at redshifts 
where well-known rest-frame optical emission lines (e.g., Hf3, 
[Olll], Ha, [Nil]) shift into the spectroscopic coverage of 
near-infrared spectrographs. 

We obtained one night from NASA on the Keck II 10m tele- 
scope with NIRSPEC to obtain near-infrared spectroscopy of 
HETDEX pilot survey LAEs. We selected LAEs from the pi- 
lot survey sample by requiring that the redshift be such that 
Ha was observable with NIRSPEC (i.e., z < 2.7) and that the 
Lya line flux was greater than 10"'* erg s~' cm~^. The lat- 
ter criterion ensured that Ha and H/3 should be bright enough 
to observe in 90 minutes with NIRSPEC (assuming Lyo/Ho 

4). Additionally, to keep the target in the NIRSPEC slit 
throughout an exposure, a star with Kvega < 18 mag is re- 
quired within ^ 38" of the target, such that both can be placed 
in the slit simultaneously. We selected ~ 10 LAEs in this 
manner, and chose the three brightest for observation on this 
run - two in the COSMOS field ( Scoville e t al. 200 7j), and 
one in the Hubble Deep Field North (HDFN; Dickinson et al.| 
|2()03), named HPS (HETDEX Pilot Survey) 194, HPS256 and 
HPS419, respectively, where the numbers correspond to the 
index in the HETDEX pilot survey catalog (Paper I). These 
three LAEs have Lya-based redshifts of 2.2889, 2.4914 and 
2.2357. More details on these LAEs are listed in Table 1. 

2.2. NIRSPEC Observations 

Our observations took place on UT 21 February 2010, and 
conditions were photometric throughout the night. We ob- 
served in low-resolution spectroscopy mode, with a slit width 
of 0.76" and the standard slit length of 42". Throughout the 
night, the cross disperser was kept at 35.76°. We observed 
with the NIRSPEC-7 filter (hereafter referred to as the K- 
band), resulting in a wavelength coverage of 2.02 - 2.43 fim, 
and the NIRSPEC-5 filter (hereafter referred to as the H- 
band), resulting in a wavelength coverage of 1.48 - 1.76 fim. 
The resolution in the K-band was R ^ 1500. 

We observed all three LAEs for 90 minutes in the K-band 
filter, obtaining 6xl5-minute exposures. We also observed 
HPS 194 for 90 minutes (6 x 15-minutes) in the H-band, while 
we only obtained 20 minutes (2x10 minutes) in the H-band 
on HPS256 due to time constraints. Telluric standards were 
observed before and after each observing setup, where we se- 
lected stars from the Hipparcos survey catalog with a spectral 
type near FOV, which minimizes hydrogen absorption lines 
common in early-type stars, and metallic lines common in 
later-type stars. Both science and calibration spectroscopic 
observations were taken in an ABBA pattern. In addition, we 
obtained arc lamp calibration images between each observ- 
ing setup, as the wavelength solution may drift throughout 
the night. NIRSPEC is known to have significant persistence 
when the integrated counts pass ^10^ per pixel. During the 
afternoon, we examined the decay of this persistence by tak- 
ing dark frames, and we found it to disappear on timescales of 
~ 15 minutes. However, as a precaution, we varied the posi- 
tion on the slit of both our objects and our standard stars, such 
that adjacent exposures should not have science data falling 
on the same detector rows. 

Our targets were acquired using the "invisible acquisition" 
mode of NIRSPEC. In this mode, an alignment star must be 



placed in the slit simultaneously with the object, and the star 
must have Kvega ^ 18 to be useful. We found alignment stars 
of Kvega = 16.5, 17.0 and 18.0 at distances of 26", 17" and 
31" from HPS 194, HPS256 and HPS419, respectively. We 
first acquired the star at the center of the slit, then we slid 
the star toward one end of the slit, such that the star and the 
target LAE were equidistant from the slit center. During an 
exposure, we manually guided the slit using the slit-viewing 
camera (SCAM), which continuously images the same field of 
view as the spectrograph, excepting the light which is trans- 
mitted down the slit. Throughout the night our seeing was 
steady at 0.6 - 0.9", comparable to the size of the slit, thus 
when the star was well centered, only the wings of the PSF 
were visible to either side of the slit. If the star started to drift 
out of the slit, we moved the slit to keep pace, using 0.5 pixel 
increments (~ 0.09"), ensuring that our object stayed as well- 
centered in the slit as possible throughout each exposure. 

2.3. Data Reduction 

We used a combination of the Keck IDL-based REDSPEC' 
package along with our own custom IDE scripts to reduce the 
data. Using REDSPEC for the initial steps, we processed each 
object/filter combination separately. For each combination, 
we used calibrations taken the closest in time, consisting of 
four calibration star spectra and two arc lamp images (one 
neon and one argon). As the flat field is not expected to vary 
during the night, we used flat fields taken during the afternoon 
for calibration. On each combination, we first ran the RED- 
SPEC task spatmap, which uses the standard star spectra 
to compute a spatial map of the image, allowing rectification 
to be performed. The task specmap was then run, which 
computes the wavelength solution, where the lines are marked 
interactively. Finally, the task redspec was run, which rec- 
tifies and extracts the spectra. However, as our objects are 
so faint that their continuum light is not visible, the redspec 
extraction was not satisfactory. Thus, we ran redspec to 
obtain the wavelength-pixel solution, but we ran the RED- 
SPEC task rectify to separately rectify the science images 
without extraction. 

We rejected cosmic rays from our rectified spectra (both 
science and calibration) using the LACOSMIC package 
(Ivan Dokkumll2001.) . The sky was then removed by subtract- 
ing adjacent image pairs, with one each in the A and B dither 
position, where the object was dithered by 5" between frames 
(with the exception of HPS419, where the dither was only 
3" in order to keep the alignment star in the slit). We then re- 
moved residual sky lines from these sky-subtracted images by 
fitting and subtracting a 15th order polynomial to each image 
column in the spatial direction, where the positive and nega- 
tive target and alignment star spectra were masked during the 
fitting. 

The ID spectral extraction was performed in IDE, by ex- 
tracting a 2D rectangular box centered on the spectra. The 
box was chosen to be 9 pixels wide, corresponding to 1.7", 
or ^ twice the average seeing. The same box size was used 
for both the target LAEs as well as the standard stars. Dur- 
ing this process, the extraction region was centered based on 
the Ha or [O III] A5007 emission lines, which were visible 
in the individual K- or H-band exposures. We also verified 
the wavelength solution during this process to find and cor- 
rect any relative shifts between the exposures (a 1.26 pixel 
shift was found and corrected in four K-band observations of 

' http://www2.keck.hawaii.edii/inst/niispec/redspec.html 



4 Finkelstein et al. 



TABLE 1 

Summary of NIRSPEC-Targeted LAEs 



Object RA Dec Vtya EWlvq Kbs (Ly") A„„t (Lya) ziya tNiRSPEC (H) tNmsPEC (K) 

(J2000) (J2000) (10"" ergs-' cm-2) (A) (A) (A) (s) (s) 



HPS194 10:00:14.18 02:14:26.11 61.0 ^ f 106+]^ 3998.0 3998.2 2.2889 ± 0.0004 5400 5400 

HPS256 10:00:28.33 02:17:58.44 31-4 ^gs 188^^^ 4244.1 4244.4 2.4914 ± 0.0004 1200 5400 

HPS419 12:36:50.13 62:14:01.07 24.4™ 72+)^ 3933.6 3933.5 2.2357 ± 0.0004 — 5400 



Note. — The equatorial coordinates provided correspond to the optical counterpart to the hya emission which were used as the NIRSPEC targets. The 
detections of rest-frame optical emission hnes at the expected wavelengths confirm these counterparts as connect (with the possible exception of HPS419). 
The corrected wavelengths are shifted to the local standai'd of rest, to correct for both the Earth's and Sun's motions at the time of the observation. The en'or 
on the redshift reflects a typical uncertainty on the central wavelength of the Lya line of 0.5 A, which dominates over the typical wavelength solution error 
of 0.05 A. The final two columns denote the NIRSPEC exposure times in each of the two filters used. 



HPS 194). These extracted regions were then summed col- 
umn by column to create a ID spectrum for each observation. 
We averaged the individual ID spectra to create the final ID 
spectrum. This procedure was repeated for each science and 
standard star observation^. 

While we were able to make ID spectra for HPS 194 and 
HPS256, we were not able to detect Ha in individual frames 
for HPS419. Given the positive detection of Ha emission in 
the other two LAEs, we would have expected a detection here. 
However, the alignment star used with HPS419 was very faint, 
with Kvega =18, thus we required 180 s SCAM integrations to 
see the star in the SCAM image. Given the level of slit shift 
seen in the other objects with much shorter SCAM exposures 
and the large distance between HPS419 and its alignment star, 
we believe that it is likely that HPS419 drifted out of the slit 
multiple times, which resulted in the non-detection. Future 
observations of very faint targets are advised to use brighter 
alignment stars (Kvega < 17.5) to avoid this issue. For the 
remainder of the paper, we focus on HPS 194 and HPS256, 
which have significant emission line detections. 

The spectra of these two LAEs were corrected for Telluric 
absorption as well as flux calibrated using the stand ard star 
spectra. This was done by taking a iKuruczl (Il993h model 
star spectrum of the same spectral type as the standard star, 
and normalizing its flux such that the H or K-band magni- 
tude of the model matched that of the standard, where the 
photometry of the standard was obtained from the Two Mi- 
cron All Sky Survey (2MASS) Point Source Catalog. We 
computed a calibration array by taking the ratio of the nor- 
malized model spectrum to the observed standard spectrum, 
interpolating over absorption features common to both spec- 
tra. This calibration array was multiplied into the observed 
object spectrum, thus performing the flux calibration and cor- 
recting for Telluric absorption. In addition, as our targets are 
unresolved from the ground, this calibration also corrects for 
point-source slit losses, as the standard and object had the 
same spectral extraction size, and the standard spectrum was 
normalized to match the total magnitude for that star 

A noise spectrum was created for each exposure for each 
object. This was done by first extracting a series of regions of 
the same size as used on the object, starting at the bottom of 
the 2D spectrum, and moving successively up by 9 rows at a 
time (i.e., so that no two rows are included in more than one 
noise estimate) in every input image, providing a ID spectrum 

- For both K-band observations and the H-band observation for HPS256, 
one standard spectrum was significantly deviant from the other three obser- 
vations; however the deviant spectrum was excluded from the analysis by the 
use of a median average. 



for each of the extracted rows. The noise in each spectral pixel 
was obtained by fitting a Gaussian to the pixel values in each 
column, performing two iterations of 3cr clipping to reject any 
real objects, and then measuring the standard deviation. This 
process created a ID error spectrum for each image. These ar- 
rays were then added in quadrature, and divided by the num- 
ber of input frames (i.e., error propagation of a mean) to create 
a final ID error spectrum for each object. Although this spec- 
trum provides a good approximation of the photometric error 
at each wavelength, it does not include uncertainties on the 
flux calibration, which are important as some of the physical 
quantities which we wish to calculate rely on ratios between 
the two filters we have observed. To include this effect, during 
the extraction of the standard star we computed a calibration 
error spectrum from the standard deviation of the standard star 
spectra in each wavelength bin. We found that the noise from 
this process was typically much less (< 1%) than the photo- 
metric noise; nonetheless this error spectrum was then added 
in quadrature to the photometric error spectrum for each ob- 
ject to include the uncertainties in the calibration process in 
our final error spectrum. We note that we have not included 
any uncertainty due to drifting of the standard star during the 
calibration exposures, as these were short and thus this error 
is likely minimal. 

3. EMISSION LINE MEASUREMENTS 

Figure[T|and Figure|2]show the reduced and flux-calibrated 
ID specti-a for HPS 194 and HPS256, respectively, with the 
ID error spectra shown in the bottom of each panel. The left 
and right panels in each figure show the region of the spectra 
around H/3 and [Olll] (left) and Ha and [Nil] A6583 (right), 
using the measured redshift of Lya to estimate the expected 
positions of these rest-frame optical lines. In Figure [T] it is 
apparent we have detected Ha emission in HPS 194 at high 
significance. We associate the weak positive feature at A = 
2.164 ^m as [Nil]. In the left panel, both [Olll] A4959 and 
[O III] A5007 are strongly detected. There also appears to be 
an emission line at the expected position of Hf3. Although 
the blue portion of this line coincides with a night sky line, 
the right portion is unaffected. One can see the red half of 
the H/3 emission line when examining the 2D spectra, thus 
we conclude it is real. In Figure |2] we again see a strong 
detection of Ha in HPS256, though there is no positive feature 
at the expected position of [N ll]. In the H-band spectrum, we 
see possible emission lines at the expected positions of [O III] 
4959, [O III] 5007 and Hf3, though it is difficult to tell if these 
lines are significant, given the increased noise present due to 
the shorter exposure time of this observation. 
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Fig. 1. — Left: H-band spectrum for HPS194. The top panel shows the reduced two-dimensional spectrum. In the middle panel, the black line shows the 
observed ID spectrum, where for display purposes we have boxcar smoothed it to match the NIRSPEC resolution at these wavelengths. Based on the redshift 
of LycK for this object, we easily identify both [O III] emission lines. We also identify H/3 emission, although the blue half of the line is affected by a sky line 
residual. The red line denotes the best fit Gaussians to these emission Unes, where the width of the Une denotes the Ict range of allowable fits. We fixed the 
redshift and FWHM of H/3 to match [O III] A5007 to allow us to fit a Gaussian to half of a line. The bottom panel shows the ID error spectrum, in the same units 
as the middle panel, magnified on the flux axis to show the structure of the noise. Right: K-band spectrum for HPS 194. While Ha is obviously strongly detected, 
[N II] only appears as a weak bump, and is not significant. Strong sky-line residuals are masked out by gray vertical bars. 
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Fig. 2. — Left: H-band spectra for HPS256, with the lines and shaded regions the same as in Figured The redshift and FWHM of the [O III] A4959 and H/3 
lines were fixed to the value measured for [O III] A5007. This observational setup only received 20 minutes of exposure time, thus this spectrum is relatively 
noisy. Right: K-band spectrum for HPS256. Although we do not see [N II] emission by eye, we still allow a line (restricted to have a non-negative line flux) to 
be fit to its expected position to obtain an estimate of the upper limit on its line flux. 



In order to measure the properties of the emission lines, 
we fit a Gaussian curve to each detected line using the MP- 
FIT IDE software package^^. For the K-band spectra for both 
objects, we fit a double Gaussian around the position of the 
detected Ha flux, allowing it to fit one Gaussian to the Ha 
flux, and a second Gaussian to [N ll]. We fit a 500 - 600 A- 
wide region of the spectra to allow ample room for the estima- 
tion of the continuum. MPFIT requires as input an estimate 
of the fitted parameters, which are the continuum flux, line 
wavelength, line full-width at half-maximum (FWHM) and 
line flux. These numbers were estimated from the ID spectra. 

As the [N II] line is either only weakly detected, or not de- 
tected at all, we forced the second Gaussian to have the same 

^ http://www.physics.wisc.edu/~craigm/idl/fitting.html 



FWHM and redshift as the first Gaussian for both objects, thus 
the only free parameter for [N ll] is the line flux. This was per- 
formed iteratively, by forcing the FWHM and redshift of the 
[N II] fit to be fixed to initially the same value as was estimated 
for Ha, and then in successive iterations forcing the [N ll] pa- 
rameters to be equal to the Ha FWHM and redshift from the 
previous iteration. This was done until the difference between 
the two Hnes was Az < 0.00001 and AFWHM < 0.01 A. In 
addition, the flux of the [N ll] line was constrained to be pos- 
itive, as this forbidden line cannot be seen in absorption. 

Similar Gaussians were fit to the H-band spectra for both 
objects. As these lines are more separated then Ha-i-[N ll] we 
fit separate Gaussians to H/3 and the two [Olll] lines. We 
fixed the ratio of the [Oinl F5nn7/F4959 to equal the theoreti- 
cal ratio of 2.98 (IStorey & ZeippenllioOOh . We fixed both the 
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TABLE 2 

Measured Emission Line Properties 



Object Line Am 


^obs 




z 


F;,„j EW,-„, 


(A) 


(A) 


(A) 




(10-" erg (A) 



HPS 194 


H/3 


4861 


15978.7 ± 0.1 


15978.2 




2.43 ± 0.43 


33 ±6 


HPS 194 


[OUT] 


4959 


16300.8 ± 0.1 


16300.3 




5.67 ± 0.48 


66 ±6 


HPS 194 


[OIII] 


5007 


16458.6 ± 0.2 


16458.1 


2.28702 ± 0.00003 


16.90 ± 0.48 


202 ±6 


HPS 194 


Ha 


6563 


21577.6 ±0.4 


21577.0 


2.28767 ± 0.00005 


8.80 ± 0.27 


177 ± 5 


HPS 194 


[Nil] 


6583 


21643.4 ±0.4 


21642.8 




< 0.15 


< 3 


HPS256 


H/3 


4861 


16971.3 ± 0.8 


16970.8 




2.24 ±0.31 




HPS256 


[Out] 


4959 


17313.5 ± 0.8 


17313.0 




<0.47 




HPS256 


[OlII] 


5007 


17481.1 ± 0.8 


17480.6 


2.49123 ±0.00017 


4.11 ±0.47 




HPS256 


Ha 


6563 


2291 1.7 ±0.44 


22911.0 


2.49094 ± 0.00007 


4.44 ± 0.29 




HPS256 


[Nil] 


6583 


22981.5 ±0.44 


22980.8 




< 0.18 





Note. — The redshifts of H/3 and [O III] A4959 were fixed to that of [O III] A5007, and the redshifts of [N II] 
were fixed to that of Ha. Both [N II] lines as well as [O III] A4959 in HPS256 were detected at < 3 cr significance, 
thus we list the 1 a upper limits. Using the redshifts listed above, the weighted mean redshifts for HPS 194 and 
HPS256 are z = 2.28712 ± 0.00002 and z = 2.49098 ± 0.00006, respectively. The equivalent widths were computed 
from the ratio of the observed line flux to that of the continuum flux near the line from the best-fit model. HPS256 
is not detected in the IRAC bands, thus the EWs for the emission lines in that object are not well constrained and so 
are not listed. 



redshift and the FWHM of the weaker H/3 and [O III] A4959 
lines to match that of the [O Hi] A5007 line during the fit- 
ting. For HPS 194, a portion of the H/3 line was lost due to 
sky emission, thus the much greater noise on the blue half 
of H/3 ensures that the red half (where we see the significant 
flux) dominates the fit. For all fits, the emission line flux was 
constrained to be non-negative. 

It is difficult to judge the detection significance of a given 
emission line simply by examining the spectra around the line, 
as surrounding regions might be heavily affected by residual 
sky lines, while the line region itself might be relatively free. 
We thus quantified the uncertainties in the emission line fit 
parameters by running a series of 10^ Monte Carlo simula- 
tions on each spectrum. In each simulation, every spectral el- 
ement was varied by a Gaussian random deviate proportional 
to the flux error at that wavelength. The altered spectrum 
then had its emission lines fit in the same way as was done 
on the actual data (i.e., holding the same parameters fixed 
which were with the data). Through these simulations, we 
thus compiled 10"' estimates of the line wavelengths (and thus 
redshifts), FWHMs and fluxes. We computed the Icr uncer- 
tainty on these parameters as one half of the spread of the 
central 68% of these values. The simulation results are shown 
in Figure [T] and Figure |2] as the width of the red line, which 
shows the central 68% of the simulation fits, highlighting the 
1(7 uncertainty on these fits. In Table 2 we tabulate all of the 
measured emission line properties, as well as the uncertainties 
on these properties from our simulations. We find that [Nn] 
is not significantly detected in either object, and [O III] A4959 
is not detected in HPS256. We detect the remaining lines at 
more than 5cr significance. 

In §7 we use upper limits on the [N ll] line fluxes to place 
constraints on the metallicities. To confirm our original es- 
timate of the 1(7 uncertainties on the [Nil] lines, we per- 
formed another set of simulations where we input into the 
K-band spectra of each object mock emission lines of vary- 
ing strengths at the position of [Nil]. We ran IVIonte Carlo 
simulations on each mock spectrum to assess the signal-to- 
noise ratio of the mock emission line detection. Plotting the 
signal-to-noise versus the mock line flux, we fit a line to the 
data, where the slope of this line is the Icr flux limit. We found 



la limits on the [Nil] flux of 1.47 x 10""^ and 1.75 x lO"'** 
erg s~' cm"^ for HPS 194 and HPS256, respectively. We use 
these limits in the following analysis. 

4. RESULTS 
4.1. Redshift and Velocity Offsets 

The simple fact that we detect Ha emission (along with H/3 
and [Olll]) near the expected wavelength confirms the origi- 
nal identification of our two objects as LAEs at z w 2.3 and 
2.5. However, by comparing the exact values of the redshift 
measured from Lya to that measured from the rest-frame op- 
tical emission lines, we can examine the kinematic state of the 
interstellar medium (ISM) in these galaxies. In many LBGs 
Lya has been found to have a slightly higher redshift than 
the systemic redshift of the galaxy by ^ 300-400 km s"', 
which coupled with measurements of blueshifted absorption 
lines implies that outflows may be pr esent in the ISMs of these 
galaxies (e.g., Shapley et al.ii2003c iSteidel et al.TOlO). 

Although these outflows appear ubiquitous in LBGs, the 
data did not exist until recently to probe for velocity offsets 
between Lya and the systemic r edshifts in the typically less 
luminous population of LAEs. iMcLinden et al.l (120101) re- 
cently measured [Olll] emission in two z ~ 3.1 LAEs. They 
found that the Lya emission lines in their galaxies were red- 
shifted by 343 ±18 and 126 ± 17 km s"' with respect to 
the [Olll] emission, which presumably hails from the sys- 
temic redshift. Although this is a small sample, if these off- 
sets are due to outflows, then the outflow velocities in one 
of these galaxies is smaller than those in LBGs. This result 
could easily be explainable, as LB Gs form stars at a higher 
rate (e.g., Papovi ch et al...2001: Shapley et al.ll2001l:lYan et"aD 
2005; Stark et al] |2009h " which could produce more intense 
outflows. However, a larger sample of LAE velocity offsets 
needs to be compiled before strong conclusions can be made. 

With our emission line redshifts, we have the data neces- 
sary to probe velocity offsets in our LAEs. However, the op- 
tical and NIR datasets were taken on different dates and at 
different observatories. The relative motion of the Earth (both 
rotational and orbital) as well as the motion of the Sun with 
respect to the observed target will thus be slightly different 
for the two datasets. Thus, before we compare the redshifts of 
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the emission lines, we need to correct the observed line wave- 
lengths for this effect. We computed the necessary correc- 
tion from the observed data to the local standard of rest using 
the online Wlsr calculator"^. The optical spectra of HPS 194 
and HPS256 were obtained at McDonald Observatory on UT 
2008-12-29 and 2008-11-29, respectively. Including the time 
of observation, we found that the radial velocity of the ob- 
server with respect to the target was -15.55 and -21.37 km 
s"', respectively. The NIR spectra of these two objects were 
obtained at the Keck Observatory on UT 2010-02-21. Simi- 
larly including the observation times, the radial velocity with 
respect to the target was 8.66 and 8.40 km s"' for the HPS 194 
H- and K-band spectra, and 9.09 and 8.94 km s"' for the 
HPS256 H- and K-band spectra, respectively. The sign of 
these measurements is such that a positive number means that 
the observer was moving away from the target (with respect to 
Vlot) at the time of the observation, so the line would appear 
at a slightly higher redshift. 

We used these velocities to correct the observed emission 
line wavelengths in both datasets, and then we computed the 
redshifts of the lines, both of which are listed in Tables 1 and 
2. For the NIR data, we only computed the redshift when it 
was a free parameter, as in both objects the redshift of [N ll] 
was fixed to match that of Ha, and in HPS256, the redshift of 
H/3 and [O III] A4959 was fixed to match that of [O III] A5007. 
We computed the systemic redshift of each galaxy using the 
weighted mean of all measured rest-frame optical emission 
line redshifts. 

The systemic redshift for HPS 194, using Ha and both 
[Olll] Unes, is z,,,, = 2.28712 ± 0.00002. For HPS256, us- 
ing Ha and [O III] A5007, it is Zsys = 2.49098 ± 0.00006. The 
Lya based redshift for these two galaxies are ziya = 2.2889 ± 
0.0004 and 2.4914 ± 0.0004, respectively. Comparing these 
two redshifts, we find that for HPS 194, Lya is redshifted com- 
pared to the systemic redshift by Av = 162 ± 37 km s"'. For 
HPS256, Lya is redshifted compared to Zsys by Av = 36 ± 
35 km s"' . However, the error on the weighted mean redshift 
does not account for systematic differences between the val- 
ues used to compute the mean. For example, in HPS 194, the 
redshifts of the Ha and [O III] A5007 lines are significantly 
different. These velocity differences are well within the in- 
strumental resolution, thus we conclude that they are a sys- 
tematic associated with the wavelength calibration. We add a 
systematic error term to the above result, where the system- 
atic error is the standard deviation of the velocity offsets from 
the individual lines used to compute the mean. Our final ve- 
locity offset results are thus: Av = 162 ± 37 (photometric) ± 
42 (systematic) km s"' for HPS 194, and Av = 36 ± 35 ± 18 
km S-' for HPS256. 

These offsets may be due to large-scale outflows in the 
ISMs of these galaxies, although only the result for HPS 194 
is statistically significant. Compared to the outflows in LBGs, 
the measured velocity offsets in our LAEs appear lesser in 
magnitude, similar to what was found in Mc Linden et al] 
( 12010) for one of their two LAEs. While a sample of four ob- 
jects is insufficient to make robust conclusions, it is intriguing 
that most LAEs appear to have smaller outflow velocities than 
LBGs. This could be due to a variety of physical effects, as we 
still do not truly understand the difference between these two 
populations of galaxies. However, in the literature, LAEs tend 
to be less luminous, lower in stellar mass, forming stars at a 

http://fuse.pha.jhu.edu/support/tools/vlsr.htmI 



lesser rate, and living in smaller dark matter halos than LBGs 
(e.g., |Papovich et al. 2001; Shaplev et al. 2003; Ga wiser et alj 
0006^ 2007; Pirzkal et al. 2007; Gronwall et al.' '20071 
Finkelstein et al. 2009d; Stark etal. 2009; Finkelstein et al] 
2010b; lOno et al.ll2010bl) . ISteidel et alj i2010) analyzed their 
sample of z ~ 2.3 galaxies to discern if the outflow veloc- 
ities they derived were correlated with any physical prop- 
erty, including those stated above. They found that the ve- 
locity differences between the systemic redshift and the inter- 
stellar absorption features and Lya emission in their sample 
do not positively correlate with any of the physical proper- 
ties. This is surprising, as it has been shown that t he out- 
flow s peed can correlate with the star-format ion rate ( Martini 
2005; Rupke et alj2005h . ISteidel et all dMol) did find that the 
velocity offset of the ISM absorption features anti-correlates 
with both the baryonic and dynamical masses (though at less 
than 3a), which could imply that the more massive galax- 
ies are accreting inflowing material at the systemic velocity. 
Thus, it appears that more work is needed to probe the root 
physical cause of the outflows in high-redshift star-forming 
galaxies. We note that internal absorption of the blue half 
of the Lya emission line can result in the observed redshift 
of Lya being s lightly hig he r than systemic witho ut invoking 
outflows (e.g.. IZheng et aH 120101 ; iLaursen et alj |2010). Fu- 
ture deep optical spectroscopy can measure the redshifts of 
absorption lines in the ISM in these LAEs, which can be used 
as a second probe for the existence of outflows. 

5. SIGNATURES OF ACTIVE GALACTIC NUCLEI 

Although star formation can produce Lya emission, active 
galactic nuclei (AGNs) also produce copious amounts of ion- 
izing photons, and thus their host galaxies emit in Lya as 
well. Understanding whether the production of Lya pho- 
tons is dominated by star formation or AGN activity is cru- 
cial in understanding not only the emission line, but also the 
characteristics of the galaxy as well. At high-redshift (z > 
3), X-ray emission and high-ionization state emission lines 
(i.e.. Civ A1549; Heii A1640; Cm] A1909) are used to 
diagnose the presence of AGNs. Using these methods, the 
AGN fractions of LA Es at z > 3 have been found to be typ- 
icaUy - 2% (e.g., Malhotra et al. 2003; Wang et al."2004j; 
iGawiser et aL.200 6; Ouchi et al. 2008) . The fraction of AGNs 
among LAEs appears to rise towards lo wer redshifts, with an 
AGN-fraction of ~ 4 - 5% at z ~ 2.1 dNilsson et al."2009t 
Guaita et al. 20 lO), and up to 40% at z ~ .3 (Finkelstein et al] 
I2009ai; iScarlataetallllOOgt ICowie et alj 12010 ). While this 
may indicate that lower-redshift LAEs are more likely to host 
AGNs, it could also indicate that less luminous, or obscured 
AGNs are lurking in LAEs at high redshift. This is especially 
critical at the redshifts of our LAE sample he re, as they he 
near to the redshift of peak AGN activity (e.g.. lOsmeij|1982l ; 
iRichar ds et al. 2006; Silverman et al. 200^. 

At lower redshifts, the mid-infrared fluxes of galaxies can 
be used to diagnose the prese nce of AGNs, as t hey typically 
exhibit red power law slopes dStern et al.ll2005l ; iDonlev et alj 
120081) . However, this typically requires detections in multiple 
mid-infrared bands, which is not always possible. Locally, 
AGNs have been identified using the line ratio diagram of 
Baldwin et al. ( 1981, hereafter BPT), which separates galax- 
ies into star-forming and AGN sequences, independent of X- 
ray luminosity. This separation is done by plotting the ratio of 
[O lll]/H/3 versus [N ll]/Ha, where AGNs experience elevated 
levels of both ionized metal lines with respect to star-forming 
dominated galaxies due to the harder UV spectra of AGNs. 



8 



Finkelstein et al. 




-2.0 -1.5 -1.0 -0.5 0.0 
Log ([Nil] A6583/Ha) 



Fig. 3. — Our two LAEs are plotted on a BPT line-ratio diagnostic dia- 
gram, with HPS 194 in red, and HPS256 in blue. The light and dark gray con- 
toured regions denote where local star-forming and AGN-dominated galax- 
ies lie, respectively, using data from the SDSS. The dotted line is the max- 
imum starburst curve from Kewley et al. (2001), while the das hed curve is 
an updated SF/AGN demarcation Hne from Kauffmann et al " f2003). The 
solid hnes denote regions expected to be inhabited by Seyfert 2's and LIN- 
ERs ( Kauffmann et al. 2003). Green triangles denote four 7, ~ 2.3 LBGs 
studied by Erb et al. 1 2006). The puiple square is the position in this dia- 
gram of the z=2.73 lensed galaxy the 8 o'clock arc ( Finkelstein et al. 20093), 
while the yellow inverted triangles represent three z ~ 2 lensed galaxies stud- 
ied bV | Hainline et al. 12009) . The gray stars denote z ~ 0.3 LAE-analogs 
IFinkelstein et aU2009bll 

To Study the possibility of AGN activity in our sample of 
LAEs, we first analy zed both the C OSMOS Chandra Bright 
Source Catalog (v2.1 iPuccetti et al. 2009) and the COSMOS 
XMM Point-like Source Catalog for X-ray counterparts near 
to our LAEs. In both catalogs (which have flux limits of ^ 5 x 
10"'^ erg s"' cm"^ in the XMM soft-band and Chandra full- 
band, respectively), we found no X-ray counterparts within 
10" of our LAEs. We then plotted the positions of our two 
LAEs on a BPT diagram, shown in Figure [3] We denote the 
limits due to the upper limit on the [N ll] fluxes by arrows. We 
plot contours o f where galaxie s from the Sloan Digital Sky 
Survey (SDSS; York et al.ll200 0) fafl in this planed with the 
AGN-dominated objects shown as the dark-gray contours, and 
the star-formation dominated objects shown as the light-gray 
contours. Both of our LAEs lie far from the AGN sequence, 
with their limits possibly pushing them even further, thus we 
conclude that the ionization in these objects is dominated by 
star-formation activity. 

Examining our two objects, HPS256 appears consistent 
with the local star-forming sequence, while HPS 194 lies 
above and t o the left, close to the maximum starburst curve of 
iKewlev et"al . (2001). A similar ti-end of enhanced [ O llll/H/3 
ratios is seen in many star-forming galaxies at z ~ 2. lErb et al] 
( 120061) studied four z ~ 2 LBGs in this plane, and found 
that even the ones far from the AGN sequence appear en- 
hanced in [Olll]/H/3 relative to the SDSS star-forming galax- 
ies (these may also be enhanced in [N II]/Hq;, moving them 
diagonally off the star-forming sequence). This effect was 
also seen in lensed LB Gs a t similar redshifts st udied by 
[Finkelstein et alj (I2009cl) and iHainUne et al] (120091) . Lastly, 
enhanced [O llll /H/3 ratios have also been noted in the z ^ 2.3 
galaxy BX418 dErb et al J 120101) which exhibits weak high- 

^ We used line fluxes from a sample of ~ 10^ SDSS star-forming galax- 
ies and ~ 3 X 10^ SDSS AGNs from the CMU-Pitt Value-Added Catalog: 
http://nvogre.phyast.pitt.edu/vac/ 



ionization emission in their rest-frame UV, similar to the 
z 2.5 galaxy MTM 095355 -545428 which also does not 
appear to host AGN activity dMalkan et al.lll996h . These 
observations together imply that the Hll regions in high- 
redshift galaxies have a stronger ionization field than seen 
locally, which could be due t o a reduced metallicity or an in- 
creased electron density (e.g., lLiu et"ani2008l:lErb et al.ll2006l: 
iBrinchmann et alj 120081) . Further data are needed to deter- 
mine if elevated [Olll]/H/3 ratios extend to the majority of 
LAEs, and to diagnose whether it is purely linked to star for- 
mation. However, we conclude that there are no indications 
that either LAE in our sample hosts an AGN. 

6. HYDROGEN LINE DIAGNOSTICS 
6.1. Dust Extinction 

We can use the spectroscopic measurements of Ha and 
H/3 emission to compute the Balmer decrement in our LAEs, 
which is a measure of how much the Ha/H/3 ratio has been 
increased due to differential dust attenuation of the bluer line. 
To perform this measureme nt, we assume an i ntrinsic ratio 
of Ha/H/? = 2.86, foUowing l Osterbrockl (119891) for the con- 
ditions in a typical H II region (and Case B re combination). 
Assuming the starburst dust extinction law of iCalzetti et alj 
(1^00), we compute the color excess due to dust extinction, 
finding E(B-V) = 0.20 ± 0.15 (Ay* = 0.8 ± 0.6 mag) and 
E(B-V) < 0.13 (Av < 0.5 mag) for HPS 194 and HPS256, re- 
spectively. The color excess for HPS256 is formally negative, 
so we quote the la upper limit here. We note that similar neg- 
ative values have been seen before (e.g., Atek et al. 2009), and 
may be indicative of H/3 photons scattering off of a reflection 
nebula within the galaxy. 

The color excess is very sensitive to changes in the Balmer 
decrement, thus the < lOcr H/3 detections result in a high un- 
certainty on the color excess, thus HPS 194 is consistent with 
Ay = as well as Ay ^ 1 mag. We note that these measure- 
ments ignore any stellar absorption, which should be a small 
effect here dRosa-Gonzalez et alj|2002h . A 1 A correction to 
the H/? EW results in an uncertainty on E(B-V) of ^ 0.05, 
much smaller than our current uncertainties. We conclude that 
higher fidelity H/3 measurements are needed to make robust 
measurements of the dust extinction levels in LAEs. Paper II 
used the UV spectral slope to estimate the level of dust extinc- 
tion in these objects, and found E(B-V) = 0.09 ± 0.05 (Ay = 
0.4 ± 0.2 mag) for HPS 194, and E(B-V) = 0.10 ± 0.10 (Ay 
= 0.4 ± 0.4 mag) for HPS256, in broad agreement with our 
Balmer decrement measurements. 

6.2. Lya Escape Fraction 

We can still gain some insight into the ISMs in LAEs by 
comparing the Ha flux to the Lya flux from these two ob- 
jects from the HETDEX pilot survey. Previous results show 
a range of Lyg/Ha of ^ 1 - 9 in LAEs at z ~ 0.3 and z ~ 
2.2 tFinkelstein et alj|2010at fHaves et al.ll2010h . We find the 
ratio of Lya flux to Ha flux in our objects to be less than 
the Case B value of 8.7, at 6.9 ± 0.5 in HPS 194, and 7.1 ± 
1.1 in HPS256^. This result immediately shows that dust if 

* Assuming a Calzetti et al. (2000) extinction curve, the conversion from 
E(B-V) to Av is a factor of ~ 4. 

' Although Lya and Ha were obtained with two different instruments, no 
"slif correction should be necessary. Lya was obtained with an IFS and thus 
is the total Lya flux, while Ha was corrected to total via the flux calibration 
(i.e., the standard star spectrum was scaled to match the 2MASS magnitude; 
see §2.3). 
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most likely present, as both ratios are significantly less than 
the Case B value, although as mentioned above, we cannot 
place strict limits on the exact value of the extinction in these 
objects. We acknowledge that the IGM can also attenuate the 
Lya flux, though significant attenuation is unlikely at z ~ 2 
(average IGM opacity at z = 2.3 would attenuate just 5% of a 
Lya line). 

If the dust extinction is very low in both objects, this implies 
that the ISM is very uniform, as only a small amount of dust is 
significantly reducing the flux ratio from the intrinsic Case B 
value. In the case of HPS 194, if the extinction is near the mea- 
sured value of Ay = 0.8 mag, this implies an inhomogeneous 
ISM, as this level of extinction would attenuate Lya signifi- 
cantly in a homogeneous geometry. The specific value of Ay 
= 0.8 mag in this object would necessitate r^ya < Tcominuum, 
implying that the Lya equivalent width (E W) we observe is 
actually enhanced over the intrinsic value ( Finkelstein et al.l 
l2009d). However, the large uncertainties on the dust extinc- 
tions in these galaxies renders us unable to make conclusions 
about the ISM, other than that a full range of geometries ap- 
pear possible. The result that both objects exhibit Lya/Ha ~ 
7 should prove useful to future surveys targeting Ha emission 
from luminous LAEs. 

Irrespective of the ISM geometry, we can estimate the frac- 
tion of Lya photons which are escaping these galaxies by 
comparing the observed Lya/Ha ratio to the Case B value. 
This would imply fe,,(Lya) = 0.80 ± 0.06 for HPS 194 and 
0.81 ±0.13 for HPS256. However, we know dust is present, 
thus these estimates are upper limits. To compute the true es- 
cape fraction of Lya photons, we must correct the Ha flux for 
dust attenuation. We do this using the values of E(B-V) from 
Paper II, which are more robust than our own Balmer decre- 
ment measurements for the reasons stated above. We assume 
that dust affects nebular emission the same as the stellar con- 
tinuum, as evidence for a preferentially higher attenuation on 
the nebular emission (Calzetti et al. 200 0& is not typically see n 
at high redshift (e.g., Erb et al. 2006; Hai nUne et all 120091) . 
though see also IFor ster Schreiberet al. (20091). We find cor- 
rected Lya escape fractions of fejc(Lya) = 0.61 ±0.11 for 
HPS 194 and fe,,(Lya) = 0.60 ± 0.22 for HPS256. These val- 
ues are consistent with those of Paper II, who find f(.jc(Lya) 
= 0.6 for HPS 194 and f„c(Lya) > 0.3 for HPS256 by 
comparing the observed Lya star-formation rates to those de- 
rived from the extinction-corrected rest- frame UV continuum. 
However, they are much higher than the median value of 
ff.sc(Lya) = 0.22 from Paper II, implying that these two LAEs 
have among the highest Lya escape fractions among galax- 
ies at z ~ 2-3. This is to be expected, as the galaxies with 
the brightest Lya fluxes were selected for NIR spectroscopic 
followup. 

6.3. Star Formation Rates 

We use our Ha observations measure the star- formation 
rates o f these two LAEs via the conversion from Kennicutt 
(119981) . which assumes a Salpeter initial mass function (IMF). 
Our observed Ha fluxes correspond to Lhq = 3.5 and 2.2 x 
10^^ erg s"', corresponding to lower limits on the observed 
SFRs of 28.0 ± 0.8 and 17.3 ± 0.9 M© yr"' for HPS 194 
and HPS256, respectively. Correcting the Ha luminosities 
for dust attenuation as in the previous sub-section, we find 
dust-corrected SFRs of 36.8 ± 5.8 for HPS 194, and 23.6 ± 
6.5 Mq yr"' for HPS256. Even without a dust correction, 
these SFRs are much larger than are typically seen in LAEs 



at z > 3 (e.g.. lGronwall et al.ll2007|). although similar to some 
bright LAEs sfiidied at z = 2.25 dNilsson et al.ll20ldh . How- 
ever, these high SFRs are again likely due to selection effects, 
as we selected the galaxies with the brightest Lya lines from 
our sample for NIR spectroscopic followup to increase the 
chances of a Ha detection, thus we also selected those with 
the highest SFRs. Additionally, as we show below, although 
our galaxies have Lya in emission, they have masses consis- 
tent with continuum-selected star-forming galaxies at these 
redshifts, which typicafly have SFRs > 10 Mq yr"' at L* 
(e.g., Reddy et al. 2006). Future work with a more homoge- 
neous sample of z ~ 2 LAEs will place stronger constraints 
on whether they are forming stars at a higher rate than LAEs 
at higher redshift. 

7. GAS-PHASE METALLICITY 

Although the presence of dust implies that many LAEs are 
not composed of pristine gas, they may be more metal poor 
than LBGs or other continuum-selected star-forming galaxies 
(i.e., BX selected) as they are less evolved in their other phys- 
ical properties. With our data, we can place constraints on 
the metallicities of high-redshift Lya-selected galaxies using 
both the N2 and 03N2 gas-phase metallicity index calibra- 
tions of (iPettini & Pagel 2004). However, [Nil] emission is 
not detected in either of our LAEs, though given the expected 
strength of the [N ll] emission line from the observed Ha flux, 
this is not surprising as an AGN-dominated ionizing spectrum 
would be required to produce a detectable amount of [N ll] 
emission. As such, we cannot measure direct metallicity val- 
ues, but we can place meaningful constraints. 

Using the 1 a upper limit on the [Nil] flux with the N2 index 
we place upper limits on the gas-phase metallicities in our 
galaxies of Z < 0.17 Zq for HPS 194, and Z < 0.28 Z© for 
HPS256 (2(7 upper limits ai-e Z < 0.25 Z© and Z < 0.41 Zq, 
respectively). The Icr upper limits from the 03N2 index are 
similar, with Z < 0.17 Zq and Z < 0.34 Zq for these two 
LAEs, respectively. 

These constraints on the metallicities of these objects are 
extremely interesting. The metallicity of HPS 194 is con- 
strained to be les s than all but one previously measured high- 
redshift galaxy. lErb et al.) (1201 Oh recently published a de- 
tailed study of the z= 2.3 BX galaxy Q2343-BX418 (hereafter 
BX418). BX418 appears unreddened, with a stellar mass of 
- 10'' Mq, and a Ha-derived SFR of 15 Mq yr"'. BX418 
also has Lya in emission, with EWrcj/ = 56 A, a Lya escape 
fraction of ^ 40%, and a somewhat broad Lya profile with 
FWHM - 850 km s"' (the Lya emission of HPS 194 has a 
similar width). They found an upper limit on its metallicity 
from the N2 index of Z^i < 0.28 Zq. They then used the 
knowledge of its low metallicity to use the lower branch of 
the R23 calibration, finding an R23 metallicity of Z«23 = 0.17 
Zq. Thus BX418 appears broadly similar to one or both of 
our LAEs in many physical properties. 

To place these metallicity limits in context, we compare 
them to the stellar masses of these galaxies, which we de- 
rive in the following section using publically available pho- 
tometry in the COSMOS field. Using the emission-line cor- 
rected SED-fitting results, we find stellar masses of 1.5 (1.0 - 
2.9) X 10'° Mq and 5.9 (2.8 - 14.2) x lO** Mq for HPS 194 
and HPS256, respectively (the values in parenthesis represent 
the 68% confidence ranges on the stellar masses). Figure]?] 
shows these two LAEs plotted on a mass-metallicity (M-Z) 
relation. We compare their positions on this plane to those 
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Fig. 4. — The positions of our two LAEs on a stellar mass - gas-phase-metallicity plane, where red represents HPS194 and blue represents HPS256. The small 
gray circles denote the SDSS sample ofTremonti et al. (2004), which represents low-redshift, star-forming galaxies, with the metallicity re-derived with the N2 
index to match the measurement used for our objects. The green triangles denote the z ~ 2.3 star-forming galaxy sample of lErb et al.l ilOOG ). which comprises 
87 LBG-like (BX-selected) galaxies split into six bins of stellar mass. The arrows for each LAE denote the 1, 2 and 3cr confidence limits on the metallicities 
of our two LAEs, while the horizontal error bar on the ia points denote the 68% confidence ranges on their stellar masses. The stellar masses are those from 
the emission-line corrected fits discussed in §8. The amount of stellar mass in old stars in HPS256 is not well-constrained, thus this object may be up to lOOx 
more massive (see §8.1). The upper limit on the metallicity of HPS194 is much lower than the z ~ 2.3 LBG M—Z trend, even though it has ~ h^y luminosity, 
implying that galaxies exhibiting Lya in emissio n may be mor e metal-poor than the general galaxy population at the same redshift. The dotted Une denotes Z = 
1/6 Zq . This is the metallicity of BX41 8 from iErb et alj j2010l) . which is the most metal-poor galaxy known at high redshift. 



of low-redshift star-forming galaxies from the SDSS taken 
from Tremonti et al. (2004), as well as LBGs at z ~ 2.3 from 
lErb et al.. (,2006) . Metallicities for both comparison samples 
have been recalculated using the N2 index. LBGs at z ~ 2.3 
lie below the low-redshift mass-metallicity sequence, though 
they follow the same trend, as lower-mass galaxies have lower 
metallicities. 

Using the metallicity limits of our objects, we study our two 
high-redshift LAEs in this plane. HPS256 is consistent with 
the z ^ 2.3 LBG M-Z relation at la. However, as we discuss 
below, the uncertainty on the stellar mass of HPS256 is high 
as it is undetected in the photometry at A > 1 fim, thus it may 
contain up to ^lOOx more stellar mass. If this is the case than 
it may reside below the z ^ 2.3 LBG M-Z relation. HPS194 
has tighter constraints on its stellar mass, and an even lower 
limit on its metallicity. We find that it lies below the z ^ 2.3 
LBG M-Z relation by at least 0.3 dex (at 2a). 

A higher ionizing field in the H II regions of high redshift 
galaxies is a possibility given the observed offsets in the BPT 
diagram of many high-redshift galaxies, including HPS 194. 
The N2 metallicity index we use was calibrated from local 
star-forming galaxies, which do not have higher-than-normal 
emission line ratios. Thus, the locally derived N2 index may 
not be directly applicable to high-redshift galaxies. The most 
direct way to show this would be to re-derive the N2 index 
with a high-redshift sample. However, this will require the 



measurement of the [O III] A4363 auroral emission line over 
a large sample, and the faint nature of this line makes this 
unfeasible with cu rrent observatorie s . How ever, the R23 in- 
dex calibration of iKobulnickv et all (119991) takes the ioniza- 
tion parameter into account, and thus may be a better alter- 
native. This is more difficult observationally, as at these red- 
shifts it will not only require J and H-band spectroscopy to de- 
tect the desired lines, it also requires a robust measurement of 
the Balmer decrement (requiring K-band spectroscopy), such 
that the emission lines can be corrected for dust extinction. 
We do not have J-band spectroscopy of our LAEs, however 
lErb et all (|2010) have computed both R23 and N2 for the z = 
2.3 BX418 galaxy, which appears to have very simil ar proper- 
ties to the LAEs in our sample. As discussed above. I&b et al.l 
(I20T0I) derive a metallicity for BX418 of Z = 0. 17 Zq with the 
R23 index, and Z < 0.28 Zq with the N2 index. The two in- 
dices produce consistent values of the metallicity, thus in this 
case at least, the offset in the BPT diagram does not appear to 
be affecting the N2-based metallicity measurements. Finally, 
even if the N2 index systematically shifts with an increased 
ionization parameter, the offset of at least one of our LAEs 
from the LBG mass-metallicity relation is real, as lErb et all 
(2006) also used the N2 index to derive the metallicities of 
their sample. 

As discussed above, SED-fitting studies have shown that 
LAEs are likely less evolved than LBGs in stellar mass, stel- 
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lar population age and dust extinction. Our above results 
show that the metallicities in at least one of our two LAEs 
are less than what one may expect from their stellar masses, 
thus LAEs may be less evolved in metallicity than LBGs, 
possibly by a large amount. This result is even more in- 
triguing given the relatively large stellar mass of HPS 194. 
From the observed photometry, this galaxy's rest-frame lu- 
minosity correspo nds to ^ 1 .4 L* us i ng the z ^ 2.3 UV lumi- 
nosity function of iReddv & Steidell (120091 the luminosity of 
HPS256 corresponds to - 0.5 L*). Thus, HPS 194 is as lumi- 
nous and massive as a typical color-selected LBG, yet it has 
a much lower metallicity. This implies that galaxies exhibit- 
ing Lya in emission may be more metal-poor than the gen- 
eral galaxy population at the same redshift, consistent with 
results found for a sample of low -redshift galaxies emitting in 
Lya (Finkelstein et al. 2010a). Although more spectroscopic 
observations of high-redshift LAEs are needed before strong 
conclusions are possible, if this result holds and LAEs are 
systematically more metal poor then the LBG population at a 
constant mass, it could confirm that LAEs are the progenitors 
of ^ L* LBGs at subsequent redshifts. 

8. EFFECT OF EMISSION LINES ON SED-FITTING RESULTS 



'Mc Linden et al.l (tlOlO*) observed that the K-band flux from 
their two z ~ 3. 1 LAEs can be wholly accounted for from their 
measured [O III] emission line fluxes, and could thus signif- 
icantly affect SED-fitting results. With our sample of LAEs, 
we can empirically examine how the presence of strong emis- 
sion lines in the rest-frame optical affects the physical prop- 
erties derived from SED fitting, using the strengths of actual 
measured lines. To perform this analysis, we used the COS- 
MOS Intermediate and Broad Band Photometry Catalog, re- 
leased in April 2009^. This catalog provides the photometry 
performed on PSF-matched images in a 3" aperture, as well 
as an aperture correction to adjust to total magnitudes. For our 
LAEs, we used the photometry from this catalog in the follow- 
ing bands (we neglected the u*-band as it is mostly blueward 
of 1216 A in both objects, and thus is subject to Une-of-sight 
IGM variations): Subaru B (0.45 /im), V (0.55 ^im), r' (0.63 
Atm), i' (0.77 ^m), z' (0.89 ^im); UKIRT J (1.25 ^im) and 
CFHT WIRCAM K, (2.15 ^m). Additionally, deep Spitzer 
Space Telescope Infrared Array Camera (IRAC) data is avail- 
able from th e S-COSMOS IRAC Photometry Catalog at 3.6 
and 4.5 ^m (ISanders et al.ll2007h . We used the S-COSMOS 
flux from a 1.9" aperture, using the provided aperture cor- 
rections to adjust to total fluxes. Although HPS256 was un- 
detected in all IRAC bands, only objects wit h > 5g detec- 
tions were included in the S-COSMOS catalog (ISanders et al.l 
120071) . Thus, for these bands we only penalized the when 
the model flux violated the 5ct limits. Additionally, while 
the catalog has a formal 2.9(7 detection in the K, band for 
HPS256, there is no significant flux apparent in the image 
(see Figure |5]l. We thus set the flux in this band to zero for 
HPS256, keeping the Icr error as the uncertainty (the stellar 
mass derived when leaving this K, flux in is consistent within 
1-2(7 of our results below). The photometry for our sample is 
listed in Table 3, and 9" cutout stamps in each band are shown 
in Figure |5] 

The photometry from these nine bands were compared to a 
suite of synthetic stellar p opulations, using the upd ated (2007 
version) of the models of iBruzual & CharloB ( 12003 , hereafter 
CB07). Models were created over a range of metallicity (0.02 

^ http://iisa.ipac.caltech.edu/data/COSMOS/datasets.htraI 



- 2.5 Zq), Star-formation history (SFH; r = 10^ - lO**'^ yr), 
age (1 Myr - iH{z\) and du st extinction (Ay = - 1.6 mag, 
using the extinction law of ICalzetti et al.l ( 120001)). IGM at - 
tenuation was applied via the prescription of iMadaul (119951) . 
The best-fit model was found via x^-minimization, with the 
stellar mass being calculated as a normalization between the 
best-fit model and the observed photometry. We added a 5% 
systematic uncertainty in quadrature to the photometric errors 
in each band to account for possible zeropoint uncertainties, 
aperture corrections, or other uncertainties which scale multi- 
plicatively with the flux. Uncertainties on the derived proper- 
ties were obtained via lO^* Monte Carlo simulations. 

The SED-fitting results are tabulated in Table 4, and the 
best-fit models are shown in Figure|6] These two LAEs appear 
to be very different physically. HPS 194 appears very old, near 
the age of the universe at z = 2.29, and is very massive, with a 
stellar mass of ~ 3 x 10'*^ Mq. HPS256 appears very young, 
with an age of only a few Myr, and is more than two orders 
of magnitude less massive. Both objects have very low stellar 
metallicity values of < 0.02 Zq, in agreement with the upper 
limits on their gas-phase metallicities. 

However, from our observed optical and NIR spectra, it is 
apparent that the emission lines we have measured in these 
objects are very strong, and they may be adversely affect- 
ing the model results. In order to quantify this effect, we 
performed a second iteration of model fitting, where we cor- 
rected the broadband photometry for the observed Lya and 
Ha emission-line fluxes. Thi s was done i n the following 
steps, following the appendix of lGuaita et al.l(l2010l) . First, we 
computed the ratio between the filter transmission at A/,„j, and 
the maximum of the transmission function, Rj-. The amount 
of flux that the line contributes to the broadband measurement 
was computed as 



fvjme -Rt ^ 



(1) 



where T(A) is the filter transmission function. This flux 
was subtracted from the catalog photometry. This correc- 
tion was made for the Z?-band for Lya accounting for a 0.30 
mag correction for HPS 194, and a 0.67 mag correction for 
HPS256. The correction due to Ha was made to the Kj-band 
for HPS 194, and accounted for a 0.18 mag correction. Ad- 
ditionally, we also excluded the J-band flux from the fit in 
HPS 194 (by making its uncertainty very large), as it appears 
to contain an excess in Figure|6l which could be due to strong 
[O II] emission, though future J-band spectroscopy is needed 
to confirm this suggestion. The correction due to the Lya 
fluxes will primarily affect the UV slope, and thus the derived 
dust extinction, while the correction due to Ha will change 
the amplitude of the z-Kg color and the derived 4000 A break, 
and thus the derived stellar population age. 

The results from the SED-fitting iteration with the corrected 
broadband fluxes are also given in Table 4, and are shown 
in Figure |6] The corrected fit for HPS256 is ^ 6 x more 
massive, and though their ages are similar, the uncertainty 
on the age in the emission-line corrected model extends to 
older ages. However, as this galaxy is not detected in its rest- 
frame optical, this fit is relatively unconstrained. Specifically, 
a single population fit cannot constrain the amount of mass in 
old stars. To examine this uncertainty, we performed a two- 
population fit, where the second population composes 96% of 
the stellar mass, and is formed in a burst at z =20. This model 
had a best-fit mass of 1.3 x lO'" Mq, which is the maximum 
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TABLE 3 

Broadband Photometry of the Sample 



Object 


B 


V 


r' 






J 




3.6 


4.5 




(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


HPS 194 


23.92 (0.06) 


24.07 (0.06) 


24.10(0.06) 


24.18 (0.07) 


23.90(0.12) 


22.73 (0.18) 


22.66 (0.14) 


22.59 (0.04) 


22.69 (0.08) 


HPS256 


24.93 (0.09) 


25.27 (0.11) 


25.39 (0.11) 


25.52 (0.14) 


25.47 (0.34) 


> 25.17 


> 25.18 


> 24.01 


> 23.32 



Note. — The optical and near-infrared data is from the COSMOS Intermediate and Broadband Photometry Catalog, while the mid-infrared data 
is from the S-COSMOS catalog. Magnitude Itr eiTors are listen in parentheses. All magnitudes are coiTected to total using the aperture connections 
from their respective catalogs. When we correct these fluxes for the measured Lya line fluxes, we find corrected magnitudes for HPS194 of B = 
24.22, and for HPS256 of B = 25.60. Likewise, correcting the measured K,-band flux for the measured Ho emission, we find K, = 22.84 in HPS 194. 
HPS256 is not detected from the J-band redward, thus we list the Icr upper limit for the J- and K,-bands, while we list the 5<j upper limit for the 
IRAC bands, as only objects brighter than these limits were included in the S-COSMOS catalog. 
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Fig. 5. — Cutout stamps (9" on a side) of HPS194 (left) and HPS256 (right) in the 9 bands that were used for SED fitting. HPS256 is undetected in all bands 
redward of 1 /im. 




Fig. 6. — The results from performing SED-fitting on HPS194 (left; in red) and HPS256 (right; in blue). The filled circles show the observed fluxes, with 
the open B and K,-band points denoting the fluxes in those bands prior to subtraction of the emission lines (an open circle also denotes the measured J-band 
flux of HPS 194, which we exclude from the emission-hne corrected fit, as it may be contaminated by [O II] emission). The hght-colored line shows the best-fit 
when emission line contamination is ignored, and the dark-colored line shows the best-fit model when the measured Lya and Ha emission fine fluxes have been 
subtracted from the B and K,-band fluxes. The squares show the bandpass-averaged flux of the best-fit emission-line corrected model. For HPS256, we show the 
Icr upper limit for the K-band flux, and the 5<t upper limits for the IRAC fluxes (as only objects brighter than this were included in the S-COSMOS catalog). In 
both objects differences arise when the emission line fluxes are removed. Specifically, in HPS 1 94 a redder UV slope is allowed after the subtraction of Ly a which 
results in lower value of tsfh and a higher extinction. Combining this with the lowered Kj-band flux results in a much younger best-fit age. The black curve in 
the right-hand panel shows the best-fit maximal mass model, which has a mass ~ lOOx greater than the single-population model, as this object is unconstrained 
at A > 1 /im. 



possible mass in this galaxy given the flux limits at longer 
wavelengths, and is a factor of > lOOx greater than the stel- 
lar mass of the single-population fit. We show this model as 
the black curve in the right-hand panel of Figure |6] 



The changes to the physical properties of HPS 194 are more 
significant. Without accounting for the emission lines, this 
object appeared to have an age of ^ 2.5 Gyr Subtracting 
the measured emission line flux from the observed photome- 
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TABLE 4 

Derived Physical Properties from SED Ftting 



Object 


Mass 


Mass 


Age 


Age 


Ay 




z 


z 


TSFR 


TSFR 




Best Fit 


68% Range 


Best Fit 


68% Range 


Best Fit 


68% Range 


Best Fit 


68% Range 


Best Fit 


68% Range 






(10"' Mq) 


(Myr) 


(Myr) 


(mag) 


(mag) 


(Zq) 


(Z©) 


(yr) 


(yr) 


HPS 194 


30.1 


28.0-31.5 


2500 


2400 - 2600 


0.08 


0.00-0.08 


0.02 


0.02 - 0.02 


109.6 


109.6 _ 109.6 


HPS 194 Corrected 


15.1 


10.3-28.7 


290 


140 - 2200 


0.24 


0.00-0.16 


0.02 


0.02 - 0.20 


10** 


106 _ io9.6 


HPS256 


0.10 


0.06-0.10 


3 


2-3 


0.00 


0.00 - 0.00 


0.02 


0.02 - 0.40 


105 


10^ - 10' 


HPS256 Corrected 


0.59 


0.28-1.42 


2 


1-57 


0.81 


0.00 - 0.97 


0.02 


0.02 - 0.02 


105 


10^ - 10' 



Note. — The best-fit physical properties from our SED-fitting analysis, along with the 68% confidence range on each property, measured from Monte 
Carlo simulations. The first row for each object shows the best-fit values for the catalog photometry, and the second row presents the results when 
correcting the photometry for the measured Lya and Ha emission line fluxes. 



try drastically reduces the best-fit age to ~ 300 Myr. Addi- 
tionally, the stellar mass is a factor of '--^ 2 x lower, and the 
derived dust extinction increases. This change can be under- 
stood when examining Figure |6] The reduction in the flux of 
the B-band results in a redder fit to the UV slope, implying a 
greater level of dust extinction''. At the same time, the less- 
ened Kj-band flux lowers the z' - color This, combined 
with the increased dust extinction, requires a much lower age 
to fit the observed z' - [3.6] color. While the difference in the 
best-fit models is substantial, the difference in mass and age 
between the two is only at the ^ 2a level, primarily because 
these objects at z ~ 2.3 are somewhat evolved, containing sig- 
nificant mass in old stars. At very high redshifts, where galax- 
ies are less evolved, and there are fewer photometric points to 
constrain the SEDs, emission line contamination likely plays 
a much larger role. 

8.1. Implications for Very High Redshifts 

The equivalent widths of our lines are difficult to measure 
from the spectra alone, as in both cases the continuum is not 
detected. The Lya EWs are computed by deriving the con- 
tinuum flux via interpolation of the flux from all detected op- 
tical filters redward of Lya, as discussed in Paper 1. For the 
rest-frame optical emission lines, we use the continuum flux 
of the best-fit emission-line corrected model to compute the 
EWs, where EW = f/,„,, / f continuum, where f continuum is the con- 
tinuum flux (in f^ units) of the model at the wavelength of 
the emission line. While the continuum flux of HPS 194 is 
well-constrained at these wavelengths due to the IRAC detec- 
tions, the same is not true for HPS256, thus only EWs for 
HPS 194 are shown in Table 2. The rest-frame EWs of both 
[O 111] lines and Ha are very large, with EW ^ 200 A for both 
[Olll] A5007 and Ha in HPS 194. 

The contributions of emission lines such as these to the in- 
tegrated fluxes of galaxies will have an even greater effect at 
high redshift, especially in the rest-frame optical. The amount 
that the emission line EWs affect the broadband photometry 
increases with redshift, as 



Am ~ -2. Slog 



1 + 



Wo(l+z) 



AA 



(2) 



due to the increase of the observed EW with 1 + z (or alter- 
natively, the shrinking of the FWHM of a filter in the rest 

' Alfliough SED fitting shows that HPS 194 appears to have Ay ~ 0.25 
mag, this is consistent with the Balmer decrement measurement, which al- 
lows a range of < Ay < 1, as well as the UV- slope derived value of Ay = 
0.4 ± 0.2 from Blanc et al. (2010). 



frame). At z ^ 7. lLabbe et al.l (12010 ) stacked a sample of 12 
z-dropout galaxies, and detect the stack in both IRAC 3.6 and 
4.5 /im bands. The resulting stellar population fit implies a 
strong 4000 A break, due to the H - 3.6 /im color The galax- 
ies used in the stack have photometric redshifts from ^ 6.5 
- 7.5, thus the [O ill] AA4959,5007 emission lines can con- 
taminate the fluxes in both of these IRAC bands. The total 
EW of both [Oni] lines measured in HPS 194 is - 270 A. 
Assuming the objects are evenly distributed in redshift, this 
will place ^ half of the line flux in each IRAC band. A EW 
^ 140 A emission line at z ^ 7 in both IRAC bands will in- 
crease the IRAC fluxes, according to Equation 1, by Ams g 
= 0.15 and Am4 5 =0.11 mag. This reduction in the ampli- 
tude of the 4000 A break of the best-fit model will result in 
a younger best-fit age. This effect is crucial, as iLabbe et alj 
620 IQ) find a best-fit age of 300 Myr, which implies that the 
stars in these galaxies began forming at z ~ 1 1 . In fact, includ- 
ing theoretical nebular emission lines in their stellar popula- 
tion models, Scha erer & de Barros (2010) find a best-fit age 
to t he photometry of Labbe et al. (2010) of only 4 Myr (see 

also lFinlatoretaLll20'lO; .Watson et al 20 IQ) . 

While measurements of objects in the mid-infrared at high 
redshift are currently very challenging, the James Webb Space 
Telescope will soon open up a wide range of high-redshift ob- 
jects for SED analyses, thus the important role of emission 
lines needs to be accounted for to ensure accurate stellar pop- 
ulation results. 

9. CONCLUSIONS 

We report detections of multiple rest-frame optical emis- 
sion lines from high-redshift Lya-selected galaxies. We de- 
tect Ha, H/3, and [O III] from two galaxies at z = 2.3 and z = 
2.5, discovered with the HETDEX pilot survey. We have used 
these emission line measurements to spectroscopically probe 
the physical properties of LAEs for the first time. 

Using the redshifts of the rest-frame optical emission lines 
as a measure of the systemic redshift, we find differences be- 
tween the redshift of Lya emission and the systemic redshift 
of Av = +162 ± 37 (photometric) ± 42 (systematic) km s"' 
for HPS194, and Av = +36 ± 35 ± 18 km s"' for HPS256. 
The slightly higher redshift of Lya is similar to what has been 
seen in more evolved high-redshift galaxies, and is potentially 
due to a large-scale outflow in the ISMs of these galaxies. 
Combining our results with the only two other LAEs which 
have had their velocity offsets measured, we find that 3/4 of 
them exhibit outflows, with most at velocities less than are 
typically seen in LBGs (> 200 km s~'). This may be due to 
lower star formation rates in the LAEs, though a larger sam- 
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pie is needed to clearly correlate the velocity offsets in LAEs 
with any physical property. 

Using the measured emission line flux ratios (and limits), 
we probe for the presence of AGN in these LAEs, using the 
BPT line-ratio diagnostic diagram. We find that both LAEs lie 
far from the AGN sequence, consistent with higher-redshift 
LAEs, in which AGNs are rare. We use the measured Ha 
and H/3 lines in our objects to measure the dust extinction. 
However, as the detection significance of H/3 in both objects 
is < 10(7, the uncertainties on the derived extinctions are high, 
thus we cannot definitively rule out zero dust in these objects. 
However, the measured ratio of Lya/Ha in both objects is ~ 
7, thus some dust is likely causing deviation from the Case 
B ratio of 8.7, though the geometry of such an ISM is as yet 
unconstrained. The escape fractions of Lya photons in these 
galaxies are high, with fejc(Lya) 60% in both LAEs. 

Lacking AGN contamination, we can use emission line ra- 
tios to measure the gas-phase metallicities. As the [N II] line 
is undetected in both objects, we can only place upper lim- 
its on their metallicities, which we find to be < 0.17 and < 
0.28 Zq for HPS 194 and HPS256, respectively. Combining 
these limits with the stellar masses measured from emission- 
line corrected SED-fitting results, we find that at least one of 
these objects lies well below the mass-metallicity sequence 
for LBGs at similar redshifts, implying that objects selected 
on the basis of their Lya emission, even if they are relatively 
massive, may be more metal poor than the general galaxy pop- 
ulation atz 2.3. 

Finally, we examine the results of SED fitting with and 
without correcting for the presence of Lya and Ha emission 
lines. We find that ignoring this strong line emission can re- 
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suit in an overestimation of the age by orders of magnitude, 
and an over-estimation of the stellar mass by a factor of two. 
This is highly significant for high-redshift studies, as the rest- 
frame EWs of both [O III] and Ha are ~ 200 A. The amount 
that an emission line affects a broadband flux is proportional 
to 1 + z, thus correcting for this contamination is crucial at 
high redshift. 

Near-infrared spectroscopy allows one to probe the physical 
characteristics of high-redshift galaxies with a much higher 
precision than is possible with SED-fitting techniques alone. 
Although it is currently time expensive, the presence of multi- 
object NIR spectrographs will alleviate this in the near future. 
The abihty of the James Webb Space Telescope to obtain mid- 
infrared spectroscopy will allow these analyses to be pushed 
to higher redshifts. 
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